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ABSTRACT 

The propagation of the solitary waves in the Bragg grating formed by 
array of thin dielectric films is considered. We assume that the thin films 
of contain the resonant molecules, which are evolved according to two-level 
atoms model, which is used to description of the coherent optical pulses 
propagation. There the alternative derivation of the Mantsyzov's equations 
is represented. 
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1 Introduction 



In during long time a special kind of gap material has been investigated £Q 
-|lUj. That was referred to as resonant Bragg grating or resonantly 

absorbing Bragg reflector (RABR) jHJ |S]. In the simplest case RABR 
consist of the linear homogeneous dielectric medium containing array of thin 
films of resonant atoms or molecules. The distance between neighboring 
films is a, the depth of film (If) is great less than wave length of radiation 
spreading in this structure. The films of two- laves interacted with ultra- 
short optical pulse have been considered and investigated in framework of 
the two-wave reduced Maxwell-Bloch equations by Mantsyzov at.al. p^- [3] , 
and by Kozhekin at.al. (HI 13 El- The existence of the 27r-pulse of self-induced 
transparency was demonstrated P3 HJ |H| • It was found [S] that the bright as 
well as the dark solitons can exist in the spectral gap, and the bright solitons 
can have arbitrary pulse area. If the two-level atoms density is very high, 
then the near-dipole-dipole interaction should be included in the Maxwell- 
Bloch equations. The effect of dipole-dipole interaction on the existence of 
gap solitons in the RABR was studied in [TO]- The many results of these 
reviewed in . Recently by using numerical simulation the optical zoomeron 
[Tlj was discovered and investigated. Optical zoomeron is similar to soliton, 
however it has the velocity which oscillates near some average value. 

The two-wave reduced Maxwell-Bloch equations are the base of the in- 
vestigation of the solitary wave propagation in the RABR. They have been 
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deduced by according to coupled-mode theory. There the alternative method 
of derivation of the evolutions equations describing the ultrashort optical 
pulses propagation in RABR was considered. Hereafter we follow the works 
PQ-[S] reasoned that the Bragg resonance corresponds for a = (A/2)m where 
m = 1, 2, 3, . . .. It will be shown that the exact equations in the form of 
discrete (recurrent) relations can be obtained. The number of approxima- 
tions after that leads to the Mantsyzov's equations. The two-wave reduced 
Maxwell-Bloch equations in more general form 0-^01 can be deduced by the 
same procedure also. 



Let us consider the ultra-short optical pulse propagation along the X-direction 
of the periodical array of thin films which are replaced in points . . . sc n _i, x n , 
x n+ i . . . (Fig.l). The medium between films characterized by the dielectric 
permittivity e. Hereafter for the sake of definiteness the TE-wave which has 
the component of the electric field parallel to the layers will be considered. 
All results can be easily generalized for the case of the TM-polarized waves. 

It is suitable the electric and magnetic strengths E , H , and polarization 
of the two-level atoms ensemble P represent in the form of Fourier integrals 



2 Transfer- operator approach 



oo 




— oo 



oo 




— oo 
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oo 



P {x n , z, t) = (27r) 2 / exp[— itot + if3z] P (x n , (3, uS)dtdz. 



—oo 



Outside off the films the Fourier components of the vectors E (x, /3, uj) and 
H (x, (3, u) are defined by the Maxwell equations. At points x n these values 
are defined from the continuity conditions. Thus, the TE-wave propagation 
can be considered in framework of the following system 



H x = -(/3/k)E, H z = -(i/k)dE/dx, E y = E, 
with boundary conditions [TTJ [T2] 

E(x n -0) = E(x n + 0) } H z (x n + Q)-H z (x n -0)=4iTrkP y (x n ,(3,u), (2) 

where k — uj/c. The solutions of equation (JTJ) in the intervals x n < x < x n+ i 
can be written as 

E(x, (3, uj) = A n (f3, uj) exp[iq(x - x n )] + B n ((3, u) exp[-iq(x - x n )], 

and 

H z (x,/3,uj) = qk~ l {A n (P,u) exp[iq(x - x n )\ - B n ((3, u) exp[-iq(x - x n ))} , 



where q = y k 2 e — (3 2 . Hence, the amplitudes A n and B n completely de- 
termine the electromagnetic field in RABR. Let us consider the point x n . 
Electric field at x = x n — 5 (5 « a) is defined by amplitudes and B^\ 
the field at x = x n + 5 is defined by A^ and B^\ Continuity conditions 
(JH) result in the following relations between these amplitudes 




(1) 




a(R) , r>(R) _ a(L) , d(I) 
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A (R) _ B (R) = A (L) _ B (L) + 4vr^ 2 g- 1 P 5 , n , 

where Ps, n — Ps{^n ^^n^) * s the surface polarization of thin film at point 
x n , which is induced by the electrical field inside film. Thus one can find 

A™ = A^ + 2 m k 2 q~ 1 P s , n , BW =B^ -2mk 2 q~ 1 P s , n . (3) 

With taking account for the strength of electric field outside off the films we 
can write 

A^l =4*)exp(iga), = B[ R > exp(-zga). (4) 

If the vectors ip ( n L) = {A£\bM) and ip ( n R) = (A^\B^) are introduced, 
then the relations can be represented in following form 

where U n is transfer operator of vector through film replaced at point 
x n . In general case it is nonlinear operator. The relations (j3J) are represented 
by vectorial form 

where linear operator of the transferring of the vector ip^f* through clearance 
between adjacent thin films V n is represented by diagonal matrix 

/ exp(iqa) 
i exp(— iqa) 

By this means we define the nonlinear transfer-operator of the vector ijj^ 
through elementary cell of RABR: 

^ L i = Wi^fntf. (5) 



In the case of linear medium transfer-operator is frequently used under con- 
sideration of the one dimensional photonic crystals, in particular, the dis- 
tributed feedback structures [T3j . 

In (JJJJ) upper index can be omitted, and the equation can be rewritten in 
form of the following recurrent relation 

A n+ i = A n exp(iqa) + 2Trik 2 q~ 1 P s , n exp(iqa) , (6) 

B n+1 = B n exp(-iqa) - 2-nik 2 q~ 1 P s , n exp(-iqa) (7) 

These recurrent relations are exact equations due to we does not use any 
approximation (for example, the slowly varying envelope of electromagnetic 
pulses approximation, the long- wave approximation). Furthermore, the sur- 
face polarization of thin film could be calculated on base of deferent suitable 
models. Here we will follow the works by B. Mantsyzov at al. PP-jS], and 
A.Kozhekin, G.Kurizki at.al. 0-05 where the two-levels atom model has 
been used. 

3 Linear response approximation 

To demonstrate that the RABR is real gap media it is suitable to obtain 
the electromagnetic wave spectrum in linear response approximation. In the 
general case we can use the following expression for polarization 

P s , n = X (u)(AW +2JW) (8) 

Substitution of this formula in ©,(0) leads to 

A n+1 = (1 + ip)A n exp(iqa) + ipB n exp(iqa), (9) 
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B n+ i = (1 — ip)B n exp(—iqa) — ipA n exp(—iqa). (10) 

Here p = p(cu) = 2irk 2 q~ 1 x(^) = 2tyujc^ 1 £ 1 ^ 2 x(^)- The wave as collective 
motion of the electrical field in grating corresponds to anzats 

A n = Aexp(ikna), B„ = B exp(ikna). (11) 

From (JHJ) , (JIOJ) it follows the linear equations are defining the amplitudes of 
the wave A and B: 

Aexp(ika) = (1 + ip)Aexp(iqa) + ipB exp(iqa), (12) 

B exp(ika) = (1 — ip)B exp(-iqa) — ipAexp(-iqa). (13) 

Nontrivial solution of this system is exists only if the determinant is not zero, 
i.e., 

/ 



det 



V 



'1 + ip) exp(iqa) — exp(ika) ipexp(iqa) 

—ip exp(-iqa) (1 — ip) exp(-iqa) — exp(ika) 



= 0. 
(14) 



Let us define 

Z = exp(ika), 

G = (1 + ip) exp(iqa) = (cosqa — psinqa) + i(pcosqa + singa) 
Then the equation (|14j) can be rewritten as equation in Z: 

Z 2 - (G + G*)Z + 1 = 0. 

The solution of that is 



Z ± = ReG±idl - (ReG) 2 . 



If ReG < 1, then ReG ± iJ 1 — (ReG) 2 = cos ka + i sin ka. Hence the 
wave numbers k± are real values and they are defined from transcendent 
equation 

cos ka = cos qa — psmqa. (15) 

That is the dispersion relation which defines the dependence of wave number 
k± versus frequency of harmonic wave propagating in linear RABR. 

If ReG > 1, then the roots of equation (fTi|l are real ones. It means that 
wave numbers k± are imaginary values. The condition Re G > 1 defines the 
frequencies of the forbidden zone. The waves with these frequencies can not 
propagate in grating. Thus the boundaries of this forbidden zone follow from 
the equation 

cosga — psinga = 1. (16) 

Model of the resonant system containing the thin films defines the explicit 
form of the function p = p(u>), that results in the different dispersion relations 
(JT5J). It should be pointed out that this dispersion relation demonstrates a 
series of gaps in the electromagnetic wave spectrum. 

4 Long- wave and weak nonlinearity ap- 
proximations 

By using long- wave approximation we can transform the exact equations ©,( 
|7|) into the differential equations. To do it is suitable to introduce the field 
variables 

A ( x ) = J2A n 5(x-x n ), B(x) = ^2B n S(x - x n ), 

n n 
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P ( x ) = Y^ P S,nK x ~ x n), 
n 

Using the integral representation for delta-function 

oo 

5(x) = (27r) _1 J exp(ikx)dk, 

— oo 

one can obtain the following expression 

oo 

A(x) = (2tv)~ 1 ^2 A n I exp[ik(x — x n )]dk = 

— oo 

oo 

= (2 I )-'Lp(ib)EA,exp(-ifa„)4. 

~, n 
— OO 

From that 

A{k) = A n exp(—ikx n ) = A n exp(-ikan) , 

n n 

and by the same way one get 

B(k) = B n exp(-ikan), P{k) = Ps,n exp(-ikan). 

n n 

The definition of these spatial Fourier components results in the period- 
icity conditions, for example, 

^(^) = y^ A n exp(-ikan) exp(±27rin) = 

n 

= y^ A n exp(—ikan ± 2nin) = (17) 

n 

= ^ A n exp [— ian{k ± 27r/a)] = ^4(fc ± 2n/a) . 

n 

From recurrent equations @,0 we have 

A(k) exp(ika) = A(k) exp(iqa) + inP(k) exp(iqa), 



B(k) exp(ika) = B{k) exp(-iqa) — ihiP(k) exp(-iqa), 

or 

[exp{ia(k - q)} - 1] A(k) = inP(k), (18) 
[exp{ia(k + q)} - 1] B(k) = -inP(k). (19) 

The linear approximation for polarization P(k) = x(co>)L4(fc) + B(k)] 
leads us to linear dispersion law once again. The system of equations 
(HEJ) , (JIH1) is equivalent of discrete equations ©,0. At this point we have 
not approximation, apart from assumption of thin films width. 

If the thin film array was absent then the dispersion relation would be 

cos ka = cos qa. 

From that, for right wave with amplitude A(k) we have k = q, whereas for 
wave with amplitude B(k) propagation in opposite direction we have k = —q. 
Let us suppose that the polarization of the thin film array produce the little 
change in wave vectors, i.e., for right wave it is k = q + 5k and for the 
opposite wave it is k — —q + 5k. Let us chouse the value of q near one of the 
Bragg resonances, to say, q = 2ir/a + 5q, where 5q -C 2Tc/a. In this case the 
equations ©j© take the form 

[exp(ia5k) - 1] A(q + 5k) = inP(q + 5k), 

[exp(ia5k) - 1] B(-q + 5k) = -inP{-q + 5k). 
By using the periodicity conditions (fTTj) these equations are rewritable as 

[exp(ia5k) - 1] A(5q + 5k) = inP(5q + 5k), 

[exp(ia5k) - 1] B(-5q + 5k) = -mP{-5q + 5k). 
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After changing of the variables 5k = 5k ± 5q we have 

~exp{ia(5k - 5q)} - l] A(5k) = mP{5k), (20) 

exp{ia(5k + 5q)} - l] B(5k) = -wP(5k). (21) 

The long-wave approximation means that nonzero values of the spatial Fourier 
amplitudes are located near zero value of argument. Hence, we can assume 
a5k -C 1 in exponential functions, that results in the following approximate 
equations 

ia(5k - 5q)A(5k) = mP(5k), (22) 

ia(5k + 5q)B(5k) = -iKP{5k). (23) 

If now we return in to spatial variable, then equations (|22|) . (J23)) lead us to 
the equations of coupled wave theory 

dA/dx = i5qA(x) + ma _1 P(x), (24) 

dB/dx = -%5qB{x) - iKa^Pix). (25) 

In these equations the both fields A(x), B(x), P(x) and parameters 5q, 
k are the functions of frequency u. To obtain final system of equations in 
the coordinate and time variables the inverse Fourier transformation would 
be done. However, frequently the assumption of a slowly varying it time 
scale envelopes of electromagnetic pulses and polarization used [Ej. This 
approximation allows simplifying the system of coupled wave equations under 
consideration. 
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5 Slowly varying envelopes approxima- 
tions 



Slowly varying envelopes approximation minds that we have to deal with 
narrow wave packets or, what is the same, with quasiharmonic waves [Ti] . 
For example, for the quasiharmonic wave the electric field is represented by 
expression 

E(x,t) — £ (x, t) exp[— iuot], 

where Uq is carrier wave frequency. The electric field E(x,t) and Fourier 
components of the envelope of the pulse E(x, t) are related by the following 
relationship 

oo 

E(x,t) = (27r) _1 J E(x,u) exp(—iut)du> = 

— oo 
oo 

= (27r) _1 J £ (x, u) exp[— i(u> + Uo)t]du> — 

— oo 
oo 

= (27r) _1 J £ (x, u — loq) exp(— icot)du>, 

— oo 

where the function E(x,u) is nonzero one if uo belongs to interval (uj — 
Alj, ljq + Au) with Au ujq. It leads to E(x, u + Uq) = £(x,u). Thus, if we 
have some relation for E(x, u), then in order to carry out the needed relation 
for £(x,u>) to be done shifting uj — > uj + u in all functions of u. 
Let be 

A(x, t) = A(x, t) exp(—iu t), B(x, t) = B(x, t) exp(—iu t), 
P(x, t) = V(x, t) exp(—iuot). 
12 



For the Fourier components A(x,u),B(x,u) and V(x,u) from (j23 [l .(|25j ) it 
follows 

dA(x, u)/dx = i5q(cu + cu)A(x, cu) + zk(uo + u)a~~ V(x, u), (26) 

dB(x, u))/dx = —i5q(ujQ + lo)B(x, uj) — in(uj + u)a~ V(x, to) (27) 

Inasmuch as in this expression A, B and V are distinguished from zero at 
uj <^ u , one can use the expansions: 

5q(cu + to) m q - 2n/a + qxuo + q 2 uJ 2 /2, k(lu + uj)a~ l w K . (28) 

where q n = d n q/duj n at u — ujq. In particular, the gf 1 = v g is the group 
velocity, q 2 takes account the group-velocity dispersion. 

With taking expansions (J28j) into account we can write following equations 
which describe the evolution of slowly varying envelopes 

< (|-^l) B+ f^ B - A * B = +jr «"- (30) 

where Ag = qo — 2Tr/a. To do next step it need to chouse the model for 
thin film medium. It can be enharmonic oscillators, the two- or three-levels 
atoms, the excitons of molecular chains, nano-particles, the quantum dots, 
and so on. Here we consider the two-level atom model. 

6 Two level atoms approximations 

The two-level atom state defines by the density matrix p. The matrix element 
p 12 describes the transition between the ground state |2) and excited state 
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|1) , p 22 and p n represents the population of theses states. Evolution of the 
two- level atom is governed by the Bloch equations [To] . 

d 

z ^^Pi2 = ^Au;p 12 - d 12 (p 22 - PnjAin, (31) 
d 

ih fo(P22 ~ Pit) = 2 ( d ^P21 A in ~ d 2 lp 12 A* n ) . (32) 

In these equations A{ n is the electric field interacting with two-level atom. 
In the problem under consideration we have Ai n = A + B and 

K oV = r - v - (p 12 ) . 

cn(u ) 

There the cornerstone brackets denote summation over all atoms within a 
frequency detuning Ao; from center of the inhomogeneously broadening line, 
ti{ujq) is the refractive index of the medium containing the array of thin 
films, n at is effective density of the resonant atoms in films. This value is 
defined thought bulk density of atoms N at , film width If and lattice spacing 
by formula n at = N at (l f /a). 

Let us suppose the group-velocity dispersion is of no importance. The 
resulting equations are the two-wave reduced Maxwell-Bloch equations. It is 
suitable to introduce the normalized variables 

ei = t di 2 A/h, e 2 = t d u B/h, x = (v g t , r = t/t . 

The normalized two-wave reduced Maxwell-Bloch equations take the follow- 
ing form 

( d d\ 

M ^7 + ^~ e i + 5e i = ~^ (P12) ' ( 33 ) 



d( dr j 
d( dr 



( d d\ 

M ^? ~ IT e 2 - 5e 2 = +7 (p 12 ) , (34) 
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d 

im Q~Pl2 = A Pl2 ~ ne in, (35) 

.9 

— n = -4Im(p 12 e*J, (36) 

where 7 = t v g /L a , 5 = t v g Aq , L a = (cn(ujo)h) / (2nuj ton a t \di 2 \ 2 ) is reso- 
nant absorption length and A = Auoto is normalized frequency detuning. We 
use the following variables n = p 22 — Pn, ej„ = e\ + e 2 . 

Let define new variables according to following expressions 

e in = e 1 + e 2 = f s exp(i5r), ei - e 2 = f a exp(iSr), p 12 = r exp(i5r). 

The system of equations (JBlIj) - can be rewritten as 

f + f = 0, (37) 

^ + ^ = 2z 7 (r) , (38) 

z^r = (A + <$)r-n/„ (39) 
d 

—n = -Alm{rf*). (40) 



From (|37jl it follows 



9/a 9/. 



<9r ' 

It allows rewrite the system ()37 p -([40 )) in another form 

d 2 f s d 2 f s n . Idr 



dC to* - 2 ^\^/' (41) 
iJU = (A + S)r-nf s , (42) 
JU = -4Im(r/;) (43) 
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If we assume that the inhomogeneous broadening is absent, i.e., the hy- 
pothesis of a sharp atomic resonant transition is true, this system is reduced 
to Sine-Gordon equation 5 + A = pQ. Injl] the steady state solution of 
(|41 )) -([43 )) with taking inhomogeneous broadening into account was found. 



7 Conclusion 

There the derivation of the equations which are describing the ultrashort 
pulses propagation in one dimensional resonant Bragg grating (or RABR) 
was considered. The grating consists in thin films array embedded into linear 
dialectical medium. The exact discrete equation for amplitudes electric field 
inside films was obtained. The long-wave approximation and slowly varying 
envelope of the pulses approximation result in the system of the coupled 
wave equations. At this point it need to chouse the model for thin film 
medium. In ^Q- [TU] the two-level atom model has been used. However, 
other models could be considered that results in new kind of RABRs. For 
example, array of thin films of the resonant optical materials with embedded 
metal nanostructures jTHj. would be described by the following system 

( d d\ 

1 [dC + d^j 61 + 661 = " 7 {a) ' (44) 

1 (lr ^) 62 ~ 662 = +1 {a) ' (45) 

d 

i—a-Aa + fi\a\ 2 a=(ei + e2). (46) 

where fi is coupling constant associated with enharmonicity of the plasmonic 
oscillations. Second example of the resonant model corresponds with three- 
level atom and two-frequency electromagnetic field. Resulting system of the 
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basic equations is similar to system considered in ^7] , where the propagation 
of the polarized waves in RABR has been studied. A fascinating model 
was proposed by A. Zabolotskii ^Hj, where thin films are replaced with J- 
aggregates. Array of thin ferroelectric films represent the gap-medium where 
the switching wave propagation could be expected. 

It should be noted that in [H] - jU] the two-wave reduced Maxwell-Bloch 
equations contain extra terms taking account the dispersion of polarization. 
These terms can be provided by this analysis also. 
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FIGURE CAPTIONS 

Fig. 1. Model of periodical nonlinear medium corresponding with resonantly 
absorbing Bragg reflector. 
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